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H ■ Abstract 

c^ ■ 

We consider Standard Model Higgs boson production through gluon-gluon fusion in 
hadron collisions. We combine the calculation of the next-to-next-to-leading order 
QCD corrections to the inclusive cross section with the resummation of multiple soft- 
gluon emissions at small transverse momenta up to next-to-next-to-leading logarithmic 
accuracy. The calculation is implemented in the numerical program HRes and allows us 
to retain the full kinematics of the Higgs boson and of its decay products. We present 
selected numerical results for the signal cross section at the LHC {^/s = 8 TeV), in the 
H — )■ 77, H — )■ WW — )■ lulu and H — )■ ZZ — )■ 4/ decay channels by using the nominal 
cuts applied in current Higgs boson searches by the ATLAS and CMS collaborations. 
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1 Introduction 

One of the main tasks of the LHC program is the search for the Higgs boson [1] and the study of 
its properties (mass, couphngs, decay widths). 

The LHC, after a successful start of pp colhsions in 2009 and 2010, has been operated at a 
centre-of-mass energy of 7 TeV in 2011, and data corresponding to an integrated luminosity of 5.7 
fb~^ have been accumulated. These data already allowed the ATLAS [2J and CMS [3] experiments 
to shrink the allowed mass range for the Standard Model (SM) Higgs boson H considerably by 
essentially excluding the Higgs bosons in the mass range (9(130 GeV) < niH < (9(600 GeV), while 
observing an excess of Higgs boson candidate events around rriH = 125 GeV. An update of the 
Tevatron results [1] with up to 10 fb~^ integrated luminosity shows a broad excess of events in the 
region 115 — 135 GeV. More data from the LHC 2012 run, that will be operated at a centre-of-mass 
energy of 8 TeV, are needed to say whether these excesses really correspond to a Higgs signal or 
are just statistical fluctuations. 

In this paper we consider the production of the SM Higgs boson by the gluon fusion mechanism 
and its decays H — )■ 77, H — )■ WW and H — )• ZZ. The gluon fusion process gg ^ H |5], through 
a heavy-quark loop, is the main production mechanism of the SM Higgs boson at hadron colliders. 
The corresponding cross section is typically at least one order of magnitude larger than the cross 
section in the other production channels (vector boson fusion, associated production....), and 
becomes comparable with the cross section for vector boson fusion only at high Higgs boson 
masses. It is thus essential to achieve reliable theoretical predictions for the gluon fusion cross 
section and the associated distributions. The dynamics of the gluon fusion mechanism is driven 
by strong interactions. Thus, accurate studies of the effect of QCD radiative corrections are 
mandatory to obtain precise theoretical predictions. 

The leading order (LO) cross section is proportional to Og, as being the QCD coupling. The 
QCD radiative corrections to the total cross section have been computed at the next-to-leading 
order (NLO) in Refs. [HI [TJ |8] and found to be of the same order as the LO contribution. The 
next-to-next-to-leading order (NNLO) corrections have been computed in Refs. P [THl [H] and 
their effect is moderate: for a light Higgs, they increases the NLO cross result by about 25% at the 
LHC {^/s = 8 TeV). We recall that all the NNLO results have been obtained by using the large- Mj 
approximation, Mf being the mass of the top quark. Corrections beyond this approximation have 
been considered in Refs. [El [131 [H [151 [IHl [IT] . 

The NNLO result mentioned above is certainly important, but it refers to a fully inclusive 
cross section. The impact of higher-order corrections generally depends on the selection cuts used 
in the experimental analysis and also the shape of the distributions is typically affected by the 
applied cuts. 

A first step in the direction of taking selection cuts into account was taken in Ref. [18], where 
the inclusive cross section with a jet veto was computed at NNLO. The first NNLO calculation of 
the Higgs production cross section that fully takes into account experimental cuts was reported in 
Ref. [19], in the case of the decay mode H — )■ 77. In Ref. [20] the calculation was extended to the 
decay mode H — > WW — )■ lulu. An independent NNLO calculation of the Higgs production cross 
section has been presented in Refs. [2T| [22] . and implemented in the parton- level Monte Carlo 
program HNNLO. Such program allows the user to evaluate the Higgs production cross section with 



arbitrary kinematical cuts and includes the decay H — )■ 77, H — )■ WW — )■ Ivlv and H -^ ZZ — )■ 4 
leptons. 

Unfortunately, fixed order calculations may suffer from perturbative instabilities when different 
energy scales are involved. An example is the transverse momentum px spectrum of the Higgs 
boson. In the small-pr region [pt ^ ^^h)i the convergence of the fixed-order expansion is spoiled 
by the presence of large logarithmic terms, dg\n^{m'\j/p'^). To obtain reliable predictions, these 
logarithmically-enhanced terms have to be systematically resummed to all perturbative orders 
[23]-[32]. It is then important to consistently match the resummed and fixed-order calculations 
at intermediate values of pt-, in order to obtain accurate QCD predictions for the entire range of 
transverse momenta. In the case of Higgs boson production the resummation has been performed 
up to next-to-next-to-leading logarithmic (NNLL) accuracy [331 EI] and matched to the fixed order 
Oi^a^ result valid at large transverse momenta [MIISHIES]- The calculation has been implemented 
in a numerical program, named HqT [37j that has been used by the experimental collaborations 
at the Tevatron and the LHC for a few years. In Ref. [38] the calculation has been extended to 
include the rapidity dependence of the Higgs boson. In Ref. [39] we have improved the calculation 
of Ref. [31] by implementing the exact form of the NNLO coefficients of Ref. |40)1 and the value 
for the coefficient A^^-* derived in Ref. [H]. The corresponding computation is implemented in a 
new version of HqT. 

In this paper we take one step forward with respect to the work of Refs. [311 [38l [39]. We 
start from the doubly differential cross section, including transverse-momentum resummation and 
rapidity dependence [38] and we implement the hard collinear coefficients of Ref. [IQ] that, together 
with the exact form of the coefficient A^'^^ [H], allow us to control the resummation at full NNLL 
accuracy. We then include the Higgs boson decay and implement the ensuing result into an 
efficient Higgs event generator, that is able to simulate the full kinematics of the Higgs boson 
and of its decay products. The resummed result is finally matched with the fixed order NNLO 
computation of Ref. [211122] to obtain a prediction that is everywhere as good as the NNLO result, 
but includes the resummation of the logarithmically enhanced contributions at small transverse 
momenta. The exact form of the second order hard-collinear coefficients of Ref. [10] permits a 
fully consistent matching with the NNLO rapidity distribution upon integration over pt- The 
calculation is implemented in a new numerical program called HRes, that embodies the features of 
HNNLO and HqT. We present a selection of numerical results that can be obtained with our program 
for Higgs boson production at the LHC (a/s = 8 TeV) up to NNLL-I-NNLO accuracy. We consider 
the decay modes H — )■ 77, H — )■ WW — )■ Ivlv and H — >■ ZZ — )■ 4/ and we compare the resummed 
results with the corresponding fixed order results, up to the NNLO accuracy, obtained with the 
HNNLO numerical code. 

The paper is organized as follows. In Sect. [2] we recall the main features of our resummation 
formalism and we introduce our NNLL-I-NNLO numerical program HRes. In Sect. |3]we present 
our numerical predictions at the LHC. In Sect. HI we summarize our results. 



'I^The results of Ref. [31] [35] were based on a reasonable approximation of these NNLO hard-collinear coefficients. 



2 Transverse momentum resummation and the HRes pro- 
gram 



We start this Section by briefly recalling the resummation formalism of Refs. [301 [311 [38] . We 
consider the inclusive hard-scattering process 



hi{pi) + h2{p2) -^ H{y,pT,mH) +X 



(1) 



where the collision of the two hadrons hi and /i2 with momenta pi and p2 produces the Higgs boson 
H with transverse momentum px and rapidity y (defined in the centre-of-mass frame) accompanied 
by an arbitrary and undetected final state X. The centre-of-mass energy of the colliding hadrons 
is denoted by y/s. 



According to the QCD factorization theorem, the doubly differential cross section for this 
process reads 



da 



dy dpr 



-{y,PT,mH,s) 



V" / dXi dX2 fai/hi{Xl,HF)fa2/h2{x2,IJ^l) 

„, n„ Jo Jo 



ai,a2 
d&ai 
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dydpj. 



{y,PT,mH,s;asifJ.R),fJ,R,fJ,l) , 



(2) 



where fa/h{x,fil) 



qf,qf,g) are the parton densities of the colliding hadrons at the factor- 



ization scale fip, d(^ab are the partonic cross sections, and /i/j is the renormalization scale. The 
rapidity, j/, and the centre-of-mass energy, s, of the partonic cross section (subprocess) are related 
to the corresponding hadronic variables y and s as: 



y = y 



1 Xi 
- In — 

2 X2 



X\X2S 



(3) 



The partonic cross section ddab is computable in QCD perturbation theory but its series ex- 
pansion in as contains the logarithmically-enhanced terms, {a^/Pr) 1ii'"(^h/Pt)) that we want 
to resum. 



To this purpose, the partonic cross section is rewritten as the sum of two terms. 



ddn 



dy dp^ 



^i.(i-es.) 

dydp^ 



+ 



da 



(fin.) 
aia2 



dydp^ 



(4) 



The logarithmically-enhanced contributions are embodied in the 'resummed' component d&a 



' a\a2 ■ 



The 'finite' component daaia2 is free of such contributions, and it can be computed by a truncation 
of the perturbative series at a given fixed order. In particular we compute daai^2 starting from 
[daaia2]io ) the usual perturbative series truncated at a given fixed order in as, and we subtract 
the perturbative truncation of the resummed component at the same order: 



"^^(fin.)" 
uu aia2 

dydp\ 


f.o. 


daa^a2 

dydp^ 


f.o. 


"^"(■•es.)" 

dydp^ 



(5) 



f.o. 



The resummation procedure of the logarithmic terms has to be carried out [23]- [2H] in the 
impact-parameter space, to correctly take into account the kinematics constraint of transverse- 
momentum conservation. The resummed component of the partonic cross section is then obtained 
by performing the inverse Fourier (Bessel) transformation with respect to the impact parameter 
b. We writqil 

-jT^{y,PT,rnH,s]as) = ^ / db - Joibpr) yVa^a2iy,b,mH, s;as) , (6) 

dydp^ s Jq 2 

where Jo{x) is the Oth-order Bessel function, and the factor W embodies the all-order dependence 
on the large logarithms ln{m'jjb'^) at large b, which correspond to ln(?7i^/p^) terms in px space. 

In the case of the pt cross section integrated over the rapidity, the resummation of the large 
logarithms is better expressed [301 EI] by defining the A^-moments Wat of W with respect to 
z = rnjj/s at fixed ttih- In the present case, in which we want to keep the dependence on the 
rapidity into account, we consider 'double' (A^i, A'^2)-iiioments with respect to the two variables 
zi = e+^mj|//Vs and Z2 = e~^m///Vs at fixed niH (note that < Zj < 1). We thus introduce 
l^(JVi,JV2) as follows M: 



^Zf'\b,mH;as)= f dz,z^^-' [ dz, z^'-' >V,,„,(y, 6,mH, S; ag) • 

Jo Jo 



(7) 



The convolution structure of the QCD factorization formula ([2]) is easily diagonalized by con- 
sidering (A?"!, A^2) -moments: 

d^iNuN.) ^ J2 fa.,h.,N.+l fa2/h2,N2-,l ^^gf ^^ , (8) 

ai,a2 

where fa/h,N = Jq dxx^~^ fa/h{,x) are the customary A^-moments of the parton distributions. 

The use of Mellin moments also simplifies the resummation structure of the logarithmic terms 
in daaia2 ^' • The perturbative factor WaiaV ^ can indeed be organized in exponential form 
as follows: 



>V(^^'^^)(6,m^;as) = T/^^^'^^^^h, «s; m^/Q^) exp{^(^-^^)(as, I>; m^/Q^)} , (9) 

where 

L^..(^ + l) . (10) 

bo = 2e~'^"^ ijE = 0.5772 ... is the Euler number) and, to simplify the notation, the dependence 
on the flavour indeces has been understood. The scale Q ~ mn in Eq. (TTOj) . named resummation 
scale, parametrizes the arbitrariness in the resummation procedure. Its role is analogous to the 
role played by the renormalization (factorization) scale in the context of the renormalization 
(factorization) procedure. Although the resummed cross section does not depend on Q when 
evaluated at all perturbative orders, its explicit dependence on Q appears after truncation of the 
resummed expression at a given logarithmic accuracy. 



■I-In the following equations, the functional dependence on the scales /i/j and hf is understood. 



The function "HC^i-^a) floes not depend on the impact parameter b and, therefore, its evaluation 
does not require resummation of large logarithmic terms. It can be expanded in powers of as as 



2 



K(^i.^^) = ao(as,mH)fl + — ?/(^-^^)«+f— ) T^^^-^^^ (^^ + . . . 1 , (11) 

L TT V TT / J 

where ao{as, mn) is the lowest-order partonic cross section for Higgs boson production. The form 
factor exp{^} includes the complete dependence on h and, in particular, it contains all the terms 
that order- by-order in as are logarithmically divergent when 6 — )■ oo. The functional dependence 
on h is expressed through the large logarithmic terms agL™ with 1 < m < 2n. More importantly, 
all the logarithmic contributions to Q with n + 2 <m <2n are vanishing. Thus, the exponent Q 
can systematically be expanded in powers of as, at fixed value of A = asi^, as follows: 

(12) 
The term Lg^^^ collects the leading logarithmic (LL) contributions oi^U^^^\ the function f^^^^ in- 
cludes the next-to-leading logarithmic (NLL) contributions as^"; g'''^^ resums the next-to-next-to- 
leading logarithmic (NNLL) terms as-^"~^, and so forth. 

Note that we use the logarithmic variable L (see Eq. (TTOl) ) to organize the resummation of the 
large logarithms ln((5^6^). In the region in which Qh 2> 1 we have L ~ ln((5^6^) and the use of 
the variable L is fully legitimate to arbitrary logarithmic accuracy. When Qb ^ 1, we have L — )■ 
and exp{^(as, iv)} — > 1. Therefore, the use of L reduces the effect produced by the resummed 
contributions in the small-6 region (i.e., at large and intermediate values of pr), where the large-6 
resummation approach is not justified. In particular, setting 6 = (which corresponds to integrate 
over the entire "pr range) we have exp{^(as, i^)} = 1: this property can be interpreted ^31j as 
a unitarity constraint on the total cross section; transverse-momentum resummation smears the 
shape of the "Pt distribution of the Higgs boson without affecting its total production rate. 

The resummation formulae ([H]), (ITT]) and flT^ can be worked out at any logarithmic accuracy 
since the functions "H and Q can be expressed (see Refs. [HI [32]) in terms of few perturbatively- 
computable coefficients denoted by A^"^ i?'^"\ iJ'^"\ CJy , G}^ , 7]y . In the case of the pr cross 
section integrated over the rapidity, Eq. Q is still valid, provided the double (A?"!, A^2)-nioments 
are replaced by the corresponding single A^-moments 'Wn-.'^n-.Qn (see Sect. 2.2 in Ref. [31] and 
Sect. 2 in Ref. [38]). 



The formalism briefly recalled in this section deflnes a systematic expansion ]3I] of Eq. Q: 
it can be used to obtain predictions that, formally, have uniform perturbative accuracy from 
the small-pT region to the large-pr region. The various orders of this expansion are denoted as 
NLL-hNLO, NNLL-hNNLO, etc., where the flrst label (NLL, NNLL, . . . ) refers to the logarithmic 
accuracy at small pr and the second label (NLO, NNLO, . . . ) refers to the customary perturbative 
order for the inclusive cross sectioi^. To be precise, the NLL-I-NLO term of Eq. ([3D is obtained 
by including the functions g^^'^ , g^^^ ^ and the coefficient H^^^ |32l [33] (see Eqs. (^ and (ITTD l 
in the resummed component, and by expanding the finite (i.e. large-pr) component up to 0{ag). 



^We note that this notation differs from the one used in Refs. [SU [351 ISSj, where the various terms of the 
expansion were denoted by NLL+LO, NNLL+NLO and so forth. Since in this paper we do not hmit ourselves 
to study the Higgs pt spectrum, we prefer to label the fixed order contributions entering the matching procedure 
according to the order they contribute to the inclusive cross section. 



(3) 



At NNLL+NNLO accuracy, the resummed component includes also the function g)^ 



^ j and 

the coefficient T-L^"^^ [40] (see Eqs. (IT2|) and ( TTTj) ). while the finite component is expanded up to 
0{ai). We point out that the NNLL+NNLO (NLL+NLO) result includes the full NNLO (NLO) 
perturbative contribution, supplemented with the resummation of the logarithmically enhanced 
terms in the small-pT region at (N)NLL. 

In order implement our calculation in a tool that can be used to perform realistic simulations, 
it is important to consider the Higgs boson decays. Since we are dealing with a scalar particle, 
the inclusion of the Higgs decay does not lead to substantial complications. However, the efficient 
generation of "Higgs events" according to the doubly-differential distribution of Eq. (jl]) and the 
inclusion of the decay are technically non trivial and require substantial improvements in the speed 
of the numerical program that evaluates the resummed cross section. The finite part in Eq. (jlj) is 
instead evaluated through an appropriate modification of the HNNLO code, which being based on 
the subtraction formalism of Ref. [H], is particularly suitable to this purpose. 



We recall [31] that, due to our actual definition of the logarithmic parameter L in Eq. Q) and 
to our matching procedure with the perturbative expansion at large pt, the integral over p^ of 
the Pt cross section exactly reproduces the customary fixed-order calculation of the total cross 
section. This integral, however, implies an extrapolation of the resummed result at large transverse 
momenta, where the resummation cannot improve the accuracy of the fixed order expansion. 
Moreover, the extrapolation of the resummed cross section at large transverse momenta may lead 
to unjustified large uncertainties and ensuing lack of predictivity (see Sect. 3 in Ref. [2I])- This is 
not a problem if the calculation is limited to the transverse momentum spectrum. In this case, in 
fact, we can simply use the fixed order result when the uncertainty of the resummed calculation 
becomes too large. In the present case, since our goal is to generate the full kinematics of the 
Higgs boson and its decays without a selection on the Higgs transverse momentum, this issue 
becomes particularly relevant. In the numerical implementation of Eq. (jlj) we thus introduce a 
smooth switching procedure at large pt, by replacing the resummed cross section in Eq. (jl]) as 
follows: 



da, 



aia2 



— )■ w{pq 



,.(res.) 
Wfaia2 



dy dpj. \ dy dpj. 

where the function w{pt) is defined aqj 



uu aia2 

dy dpip 



+ {1-w{pt)) 



da, 



aia2 



dydp\, 



(13) 



f.o. 



w{pt) 



1 Pt < p^'t' - ^Pt 
fiPr) Pt' - ^Pt < Pt < Pt' + ^Pt 
pt> p't' + ^Pt 



(14) 



and the function /{pt) is chosen in such a way that w{pt) and w'{pt) are continuous in all the 
range of transverse momenta. In particular, we choose 



f^- 



fiPTl = - COS vr 



Apj 



+ 1 



(15) 



We have checked that the parameters p^' and ApT can be consistently chosen so as not to spoil 
our unitarity constraint, and that the integral of our NLL+NLO and NNLL+NNLO resummed 
result still reproduces well the NLO and NNLO inclusive cross sections (see Sec. 3). 



^We note that a simpler switching option is available in the new version of HqT 



We have implemented our calculation in a numerical program called HRes, by considering three 
decay modes: H — )■ 77, H — )■ WW — )■ Ivlv and H — )■ ZZ — )■ 4 leptons. In the latter case the 
user can choose between H — )■ ZZ — )■ fi'^fi~e~^e~ and H — )■ ZZ — )■ e~^e~e~^e~, which includes the 
appropriate interference contribution. The program can be downloaded from [37j, together with 
some accompanying notes. 



3 Results 



3.1 Preliminaries 

We consider Higgs boson production at the LHC (e.g. pp collisions at ^/s = 8 TeV). In order to 
avoid a multiple presentation of similar results we use MSTW2008 parton distributions [H], with 
densities and as evaluated at each corresponding order, i.e. we use (n + l)-loop as at N^LL+N^LO 
and N^LO (with n = 1,2), and 1-loop as for LO. Unless stated otherwise, renormalization, 
factorization and resummation scales are set to their default values, hr = hf = 2Q = niH- We 
remind the reader that the calculation is performed in the Mt — )■ 00 limit. 

As for the electroweak couplings, we use the scheme where the input parameters are Gp, 
mzi my/ and a(m^). In particular we take Gp = 1.16639 x 10~^ GeV~^, niz = 91.188 GeV, 
mw = 80.419 GeV and a{mz) = 1/128.89. The decay matrix elements are implemented at Born 
level, i.e., radiative corrections are completely neglected. The Higgs boson is treated in the narrow- 
width approximation, but in the W and Z decays we take into account finite width effects, by 
using Tw = 2.06 GeV and Tz = 2.49 GeV. As explained in Sect. [2], in order to obtain meaningful 
predictions in the entire range of transverse momenta, we apply a smooth switching procedure (see 
Eq. (IT3II ). In our numerical implementation the parameters in Eq. (ITSl) are phenomenologically 
chosen to be Apr = 30 GeV and p?^ = amn + h{mH/Q + c)\/s. At NLL+NLO accuracy we set 
a = 1/2, 6 = 1.2 X 10-3 and c = 2.5, whereas at NNLL+NNLO we set a = 0.6, 6 = 2.8 x lO'^ 
and c = 0. We postpone some comments on the dependence of our results on these parameters to 
Sect. EH 



3.2 // ^ 77 

We first consider the production of a SM Higgs boson with mass rriH = 125 GeV. The width is 
computed with the program HDECAY [15] to be Th = 4.15 MeV. With this choice of rriH, the 
effects of finite width can safely be neglected. 

When no cuts are applied, the px spectrum of the Higgs boson obtained with HRes must 
be in agreement with the one obtained with the HqT numerical program. In Fig. [T] we compare 
the two spectra to check that this is indeed the case, within the statistical uncertainties. The 
corresponding inclusive cross sections are reported in Table [H where we show the new resummed 
results obtained through the HRes code, and we compare them with the fixed order predictions 
obtained with the HNNLO code. We see that the NLL+NLO (NNLL+NNLO) inclusive cross section 
agrees with the NLO (NNLO) result to better than 1%. 



1.000 



0.500 -- 



> 

O 0.200 



^0.100 

ft 
■a 
"\ 

■o 0.050 



0.020 



0.010 



0.005 



^ 1 1 1 


' 


1 


III 


1 1 


III 


1 


-- "% 










- HqT2.0 


— 




- 


\ 






I 


HRes 


- 


— 






sw 






— 


pp^H+X^77+X 
Vs=8 TeV, MSTW2008 


^ 


^^ 


i'^S^C^ 


— 


— Mr=MR=2Q 


=mH= 


= 125 


GeV 




^^ 




1 1 1 1 


1 


1 


1 1 1 1 


1 1 


ill! 


1 -■ 







50 



100 
Pt (GeV) 



150 



200 



Figure 1: Transverse momentum spectrum for the iJ — > 77 signal at the LHCforrriH = 125 GeV, 
obtained at NNLL+NNLO with HRes compared to the corresponding result from HqT. The result 
from HqT is multiplied by the branching ratio BR{H — )■ 77) = 2.245 x 10^^ l4^ - 

As an example, we apply the following cuts on the photons. For each event, we classify the 
photon transverse momenta according to their minimum and maximum value, PTmm and PTmax • 
The photons are required to be in the central rapidity region, |?7| < 2.5, with prmin > 25 GeV and 
PTmax > 40 GeV. Note that an isolation cut on the photons is generally required. For example, a 
standard isolation is to require the total transverse energy in a cone of a given radius R around 
each photon to be smaller than a fraction of the photon pt- Such cuts cannot be taken into account 
in our resummed calculation, since we are inclusive over the QCD radiation recoiling against the 
Higgs boson. Their effect can be estimated with the HNNLO code and turns out to be rather small. 



Cross section 


NLO 


NLL+NLO 


NNLO 


NNLL+NNLO 


Total fb 


30.65 ± 0.01 


30.79 ± 0.03 


38.47 ± 0.15 


38.41 ± 0.06 


With cuts [fb 


21.53 ± 0.02 


21.55 ± 0.01 


27.08 ± 0.08 


26.96 ± 0.04 


Efficiency [%] 


70.2 


70.0 


70.4 


70.2 



Table 1: Fixed order and resummed cross sections forpp — )■ H 
and after geometrical acceptance cuts. 



X -> 77 + X at the LHC, before 



We recall that the resummation does not affect the total cross section for the Higgs boson 
production, but when geometrical cuts are applied, their effect can act in a different way on fixed 
order and resummed calculations. In Table [1] we compare the accepted cross sections, obtained 
by the fixed order and resummed calculations, and the corresponding efficiencies. The numerical 
errors estimate the statistical uncertainty of the Monte Carlo integration. Comparing resummed 
and fixed order predictions, we see that there are no substantial differences on the accepted cross 



section, due to the fact that the integration is performed over a wide kinematical range. In 
Table [2] we report the accepted cross section for different choices of the scales. After selection 
cuts, the scale uncertainty is about ±15% (±18%) at NLL+NLO (NLO) and ±9% (±10%) at 

NNLL+NNLO (NNLO). 



Cross section fb 


NLO 


NLL+NLO 


NNLO 


NNLL+NNLO 


{2Q = fiF = I^r) = mH/2 


25.92 ± 0.02 


25.57 ± 0.03 


29.52 ± 0.13 


29.59 ± 0.11 


{2Q = fiF = m) = ITT-H 


21.53 ± 0.02 


21.55 ± 0.01 


27.08 ± 0.08 


26.96 ± 0.04 


(2Q = iiF = I^r) = 2mH 


18.17 ± 0.01 


18.80 ± 0.02 


24.43 ± 0.06 


24.69 ± 0.06 



Table 2: NLO, NLL+NLO, NNLO and NNLL+NNLO accepted cross sections for pp ^ H + X -^ 
77 + X at the LHC, for different choices of the scales. 

In Fig. [2] we study the distribution in the azimuthal separation of the photons in the transverse 
plane, A(f). At LO the photons are back-to-back, and thus A0 is 180°. Beyond LO, events with A(j) 
different from 180° are allowed, but NLO (dots) and NNLO (dashes) results show an unphysical 
behaviour as A0 — > 180°. The resummed NLL+NLO and NNLL+NNLO results lead instead to 
a smooth behaviour in this region. On the other hand, A0 — )■ corresponds to a kinematical 
configuration where the diphoton system is produced with large transverse momentum, so the 
result is fully dominated by the corresponding fixed order calculation. 
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Figure 2: A0 distribution from the ii/^ — )■ 77 signal at the LHC, obtained by the fixed order and 
resummed calculation. 

An interesting observable, which has been used by ATLAS to split the ff — )■ 77 analysis in 
categories |1S], is the thrust transverse momentum prUII l47j. Defining the thrust axis i and the 



"in the context of Drell-Yan lepton pair production, this variable is also called ar [47ll48l 



transverse momentum of the diphoton system •p^'^ as follows 

i 



Pt 



Pi' 



\pi'-Pi'y 

the pTt is then calculated according to: 

PTt 



Pi' 



pi' +Pt, 



Wi"' X i\. 



(16) 



(17) 



In Fig. Owe report the pTt distribution, obtained at NLO (dots), NNLO (dashes), NLL+NLO 
(dot dashes) and NNLL+NNLO (solid). We see that in the high pxt region the NLL+NLO 
prediction agrees with the NLO one, and the NNLL+NNLO prediction agrees with NNLO. In 
the low PTt region the NLO result diverges to +00, whereas the NNLO diverges to —00. Such 
behaviour is analogous to the behaviour of the pt distribution of the Higgs boson when computed 
at fixed order in QCD perturbation theory. The NLL+NLO and NNLL+NNLO results obtained 
with HRes are instead finite as pxt -^ 0, approaching a constant value. 
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Figure 3: pxt distribution for the iJ — )■ 77 signal at the LHC, obtained at NLL+NLO and 
NNLL+NNLO compared to the corresponding NLO and NNLO results. 



In Fig. 4(a)|4(c) we plot the photon pt distributions prmm and prmax- These distribution are 
enhanced when going from LO to NLO to NNLO according to the increase of the total cross section. 
We note that, as pointed out in Ref. [21], the shape of these distributions sizeable differs when 
going from LO to NLO and to NNLO. In particular, at the LO the two photons are emitted with 
the same pt because the Higgs boson is produced with zero transverse momentum, hence the LO 
PTmin and pTmax are exactly identical. Furthermore the LO distribution has a kinematical boundary 
at Pt = mn/^ (Jacobian peak), which is due to the use of the narrow width approximation. Such 
condition is released once extra radiation is accounted for. Thus higher order predictions suffer of 
perturbative instabilities, i.e. each higher-order perturbative contribution produces (integrable) 
logarithmic singularities in the vicinity of that boundary, as explained in Ref. [19] . 



The same prmm and PTmax predictions are shown in Fig. 4(b)|4(d) in this case the NNLO 



result is compared with the resummed result at the NLL+NLO and NNLL+NNLO accuracy. 
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Figure 4: Distributions in prmin (a,h) and pTmax (c,d) for the i/ — )■ 77 signal at the LHC, obtained 
by fixed order (a,c) and resummed (b,d) calculations. In the right panels the fixed order NNLO 
result is also shown for comparison. 



As expected [19], resummed results do not suffer of such instabilities in the vicinity of the LO 
kinematical boundary; the resummed distributions are smooth and the shape is rather stable when 
going from NLL+NLO to NNLL+NNLO0. 

An analogous perturbative instability is present in the PTmin distribution around PTmin ~ 40 
GeV at NLO and NNLO (see Fig. 4(a)). Such instability, which is not related to the use of the 
narrow width approximation, is due to the choice of asymmetric cuts for the photons. Beyond 
LO, the region pTmin < 40 GeV opens up, and the step-like behaviour at LO leads to integrable 
logarithmic singularities at NLO and beyond. The resummed NLL+NLO and NNLL+NNLO 



**We note that the smah-pT resummation we perform in this paper is not strictly the one needed to cure these 
logarithmic singularities [l^. Nonetheless, since our resummation provides a correct description of the Higgs boson 
kinematics, we do not expect that a rigorous treatment of these logarithmic singularities would lead to substantial 
numerical differences. 
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results are free of such perturbative instability. 

Finally, a variable that is often studied is cos 9* , where 6* is the polar angle of one of the photons 
in the Higgs boson rest frame. Given the 4-nionientuni of the photon p^ = {mH/2,f)^ 
Higgs rest frame, the 6* angle is defined as follows 



^ ,p1) i^ th^ 



COS0* 



bll 



mH/2' 



1^ 



considering the on-shell condition for the photon pJp +p2'^ = {mH/'^Y ^^'^ i^ak at the LO the px 
of the Higgs boson is zero, we can invert the on-shell condition, obtaining 



cosri 



Apl' 



m 



(19) 



H 



A cut on the photon transverse momentum pj implies a maximum value for cosO* at LO. For 
example for mn = 125 GeV and p^ > 40 GeV we obtain 



pI > 40 GeV 



cosri < Icos^: 



0.768. 



(20) 



At the NLO and NNLO the Higgs transverse momentum is non vanishing and events with 
I cos 0*1 > I cos0*^t| are kinematically allowed. In the region of the kinematical boundary higher- 
order perturbative distributions suffer of logarithmic singularities (as it happen for the photon 
distributions discussed above). In Fig. [5] we report both the distributions (normalized to unity) 
obtained by fixed order and the resummed calculations. We see that the resummed results are 
smooth in the region around the kinematical boundary. Away from such region, fixed order and 
resummed results show perfect agreement. 
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Figure 5: Normalised cos6* distribution at the LHC. On the left: LO, NLO and NNLO results. 
On the right: resummed predictions at NLL+NLO and NNLL+NNLO accuracy are compared with 
the NNLO result. 
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3.3 H -^ WW -^ lulu 



We now consider the production of a SM Higgs boson with mass mn = 140 GeV. The width is 
computed with the program HDECAY [15] to be T^ = 8.11 MeV. We consider the decay W ^ Iv 
by assuming only one final state lepton combination. In order to isolate the possible signal some 
acceptance cuts are needed. Here we apply the following set of cuts [50] : 



• the event should contain two opposite charged leptons having px larger than 20 GeV and in 
the central rapidity region |?7| < 2.5; 

• the missing px of the event should be larger than 30 GeV; 

• the invariant mass of charged leptons should be larger than 12 GeV. 

In Table [3] we report the fixed order and resummed predictions for the total and accepted 
cross section. We see that, also in this case, the inclusion of resummation does not lead to 
substantial differences on the accepted cross section. In Table H] we report the accepted cross 
section for different scales. The scale uncertainty is about ±17% at NLL+NLO and NLO, whereas 
at NNLL+NNLO and NNLO it is reduced to ±10%. 



Cross section 


NLO 


NLL±NLO 


NNLO 


NNLL±NNLO 


Total fb 


61.58 ± 0.04 


61.58 ± 0.04 


76.94 ± 0.09 


76.88 ± 0.19 


With cuts fb 


20.98 ± 0.03 


20.90 ± 0.02 


26.44 ± 0.10 


26.32 ± 0.05 


Efficiency % 


34.0 


33.9 


34.4 


34.2 



Table 3: Fixed order and resummed cross sections for pp -^ H + X -^ WW ± X — )■ lulu ± X 
before and after selection cuts. 



Cross section [fb] 


NLO 


NLL±NLO 


NNLO 


NNLL±NNLO 


(2(5 = iiF = /i/j) = rnH/2 


25.14 ± 0.03 


25.13 ± 0.04 


29.16 ±0.17 


29.05 ± 0.22 


{2Q = /iF = I^r) = mn 


20.98 ± 0.03 


20.90 ± 0.02 


26.44 ± 0.10 


26.32 ± 0.05 


(2Q = ^iF = Hr) = 2mH 


17.76 ± 0.02 


18.26 ± 0.03 


23.85 ± 0.07 


24.14 ± 0.10 



Table 4: Fixed order and resummed accepted cross sections for pp -^ H+X — )■ WW+X — > lulu+X 
at the LHC, for different choices of the scales. 

For each event, we classify the transverse momenta of the charged leptons according to their 
minimum and maximum value (as we did for the photon transverse momenta in if — )■ 77). 
In Fig. [6] we plot the corresponding distributions. We compare the resummed NLL±NLO and 
NNLL±NNLO predictions with the corresponding fixed order predictions at the LO, NLO and 
NNLO accuracy. We see that QCD corrections tend to make the distributions harder. Analogous 
effects are observed on the average transverse momentum spectrum of the W bosons, which is 
reported in Fig. [71 In particular, in order to quantitatively estimate the impact of the resummation. 
Figs. O [7] are organised in two panels. In the upper panels, we show the predictions obtained by 
different fixed order and resummed calculations. In the lower panel we plot (in red) the ratio 
NLL±NLO/NLO and (in blue) NNLL±NNLO/NNLO. From Fig. E we note that, in the peak 
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region, for both the PTmin and PTmax distributions, the resummed result is smaller by 5 — 10% 
at NLL+NLO and by 2 — 4% at NNLL+NNLO with respect to the corresponding fixed order 
prediction. In the intermediate region the resummation affects the results in the opposite direction, 
enhancing the cross section up to about 30% at NLL+NLO and 10% at NNLL+NNLO. The effects 
observed for the average pt of the W bosons (see Fig. [Tj) are even more pronounced. These effects 
on the Pt spectra imply that the agreement between resummed and fixed order predictions we 
have observed in Table [3] cannot persist in general. When more restrictive cuts on the transverse 
momenta are applied, we anticipate non negligible effects from resummation. 



Vs=8 TeV. MSTW2008 
MF=MH=2Q=mH=140 GeV 




Figure 6: Transverse momentum spectra of the lepton with minimum (left) and maximum (right) 
Prforpp-^H + X-^ WW + X -)■ Mu + X at the LHC. Resummed results at NLL+NLO and 
NNLL+NNLO accuracy are compared with fixed order predictions at LO, NLO and NNLO. The 
lower panels show the NNLL+NNLO result normalized to NNLO (solid) and the NLL+NLO result 
normalized to NLO (dashes). 



A very important discriminating variable for the H — )> WW — )■ lulu decay channel is the 
azimuthal separation of the charged leptons in the transverse plane, A0. As is well known |51j . 
for the Higgs boson signal the leptons tend to be close in angle, thus the bulk of the events is 
produced at small A0. Our results for the A(j) distribution are reported in Fig. |H1 We can see that 
in the very small A0 region (A0<3O°), there are less events than expected: this is an effect of 
the applied cuts. We notice that the steepness of the distributions increases when going from LO 
to NLO and from NLO to NNLO, and also increases when going from fixed order to resummed 
predictions, i.e. from NLO to NLL+NLO and from NNLO to NNLL+NNLO. This fact can be 
interpreted as follows: when the Higgs boson px distribution is harder the final state leptons tend 
to be more boosted in the transverse plane and thus their transverse angular separation becomes 
smaller. As a consequence the steepness of the A0 distribution increases and the efficiency of cuts 
slightly increases with the perturbative accuracy. 



3.4 H ^ ZZ 



IT n e^e 



We now consider the production of a Higgs boson with mass mu = 150 GeV. The width is 
computed with the program HDECAY [45j to be Th = 16.9 MeV. In this mass region, the 
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Figure 7: Average transverse momentum spectrum of the W bosons for pp — )■ H+X — )■ WW+X — )■ 
lulu + X at the LHC when cuts are applied. Resummed results at NLL+NLO and NNLL+NNLO 
accuracy are compared with fixed order predictions at LO, NLO and NNLO. The lower panel 
shows the NNLL+NNLO result normalized to NNLO (solid) and the NLL+NLO result normalized 
to NLO (dashes). 
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Figure 8: Same as in Fig. [^ hut for the A(f) distribution. 



15 



H -^ ZZ -^ 4/ decay mode is not the dominant decay channel, but still it can provide a clean 
and useful four lepton signature. In the following we consider the decay of the Higgs boson in two 
different lepton pairs. 

We consider the following cuts [50] : 

• the event should contain two pairs of opposite charged leptons 

• each lepton must have pt larger than 5 GeV and should be in the central rapidity region 
\7]\ <2.5; 

• for each lepton pair, the closest {mi) and the next-to-closest (7712) to mz are found; then 
mi, 7712 are required to be 7ni > 50 GeV and 7712 > 12 GeV. 

Note that an isolation cut on the leptons is generally required. For example, a typical isolation is 
to require the total transverse energy Et in a cone of a given radius R around each lepton to be 
smaller than a fraction of the lepton pt- As in the i7 ^ 77 decay mode isolation cuts cannot be 
applied, because in the resummed calculation we are necessarily inclusive over the QCD radiation 
accompanying the Higgs boson. By using the fixed order HNNLO code we have checked that the 
numerical effect of the isolation cuts is extremely small. 

In Tab. [5] we compare the effects of cuts on the inclusive cross sections. As in the if — )■ 77 
and H — > WW decays, the efficiency slightly improves increasing the perturbative accuracy, but 
no substantial effects from resummation are observed. In Tab. [6] the accepted cross section for 
different choices of the scales is reported. At the NLL+NLO (NLO) accuracy the ensuing scale 
uncertainty is ~ ±15% (±17%) and at the NNLL±NNLO (NNLO) it is ~ ±9% (±10%). 



Cross section 


NLO 


NLL+NLO 


NNLO 


NNLL+NNLO 


Total [fb] 


1.720 ± 0.001 


1.720 ± 0.002 


2.142 ± 0.004 


2.156 ± 0.006 


With cuts [fb] 


1.127 ± 0.001 


1.136 ± 0.001 


1.413 ± 0.005 


1.427 ± 0.003 


Efficiency [%] 


65.6 


66.0 


66.0 


66.2 



Table 5: Fixed order and resummed cross section for pp ^ H 
cross section before and after geometrical acceptance cuts. 



X ^ZZ + X ^ fi+l^-e+e- + X 



Cross section [fb] 


NLO 


NLL+NLO 


NNLO 


NNLL+NNLO 


{2Q = I2F = ij,r) = mH/2 


1.350 ± 0.001 


1.350 ± 0.004 


1.572 ± 0.007 


1.570 ± 0.006 


{2Q = fiF = JJ-r) = rriH 


1.127 ± 0.001 


1.136 ± 0.001 


1.413 ± 0.005 


1.427 ± 0.003 


{2Q = I2f = ij,r) = 2mH 


0.954 ± 0.001 


0.992 ± 0.003 


1.273 ± 0.003 


1.310 ± 0.003 



Table 6: Fixed order and resummed accepted cross sections for pp -^ H + X ^ ZZ + X — )■ 
li~^fi~e~^e~ + X at the LHC, for different choices of the scales. 

In Fig. |9]we plot the four pt spectra of the final state leptons. Note that at LO the Pti,Pt2 
are kinematically bounded by mH/2, whereas prs < mu/'i and pt4 < '"^///4. In the vicinity of 
such boundaries, higher order QCD predictions may in principle develop perturbative instabilities. 
On the other hand, contrary to what happens in the if — )■ 77 decay mode, the LO distributions 
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smoothly reach their kinematical boundary and we do not observe perturbative instabihties beyond 
the LO. The impact of resummation is to make the transverse momentum spectra harder. The 
resummation effects are more pronounced in the leading lepton transverse momentum spectrum 



(see Fig. 9(a)) and less evident in the softest lepton spectrum (see Fig. 9(d)). 
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Figure 9: Transverse momentum spectra of the final state leptons for pp ^ H + X ^ ZZ + 
X — )■ fi^fi~e^e~ + X at the LHC, when cuts are applied. The lepton pj- are ordered according 
to decreasing pt- They are obtained through fixed order (black) and resummed (red and blue) 
calculations. The lower panels show the ratios between resummed and fixed order predictions. 

In Fig. [in] we show the average pt distribution of the two Z bosons. The comments are 
analogous to those for previous distributions: QCD radiation tends to make the distribution 
harder and the fixed order results are again recovered at large transverse momentum. 
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Figure 10: Average px spectrum of the Z bosons for pp —^ H + X -^ ZZ + X — ;■ fi^fi^e^e^ + X 
at the LHC, when cuts are applied. Resummed results at NLL+NLO and NNLL+NNLO accuracy 
are compared with fixed order predictions at LO, NLO and NNLO. The lower panel shows the 
NNLL+NNLO result normalized to NNLO (solid) and the NLL+NLO result normalized to NLO 
(dashes). 

3.5 Discussion 



As explained at the end of Sec. 2 our numerical program implements a smooth switching procedure 
between the resummed and fixed order results (see Eqs. (IT3|) and (TT5|) ). The numerical parameters 
in Eq. (fT3|) can be consistently chosen such that the integral of our NLL+NLO and NNLL+NNLO 
resummed result still reproduces well the NLO and NNLO inclusive cross sections. We have studied 
the uncertainty inherent in such switching procedure. We find that for other forms of the switching 
function and different choices of the corresponding parameters the numerical effects on the final 
results remain within the uncertainties estimated by performing variations of the resummation 
scale Q. 

We conclude this Section by adding few comments on the work of Ref. [52]. In this paper the 
RESBOS generator [53], which is based on the classical 6-space resummation formalism of Ref. |27j . 
is used to perform a study of transverse momentum resummation effects in the H — )■ WW — )■ lulu 
and H — )■ ZZ — )■ 4/ channels at the Tevatron and the LHC. The resummed calculation in the low 
Pt region is matched to the 0{a^) result at high pt- Besides the differences in the resummation 
formalism (see Ref. [31] for a detailed discussion) there are a few differences with respect to the 
work presented here. Our calculation implements the value of the coefficient A^^^ from Ref. [H], 
whereas in Ref. [S2] the authors use the result of Ref. [M] that applies to threshold resummation. 
The calculation of Ref. [52] does not include the hard coUinear coefficients Ti^'^^ presented in 
Ref. [10] and thus its accuracy, with our notations, is essentially limited to NLL+NLO (plus some 
of the NNLL terms). Finally, the calculation of Ref. [52] does not exploit a unitarity constraint 



on the total cross section, and thus the normahzation of the ensuing resummed spectra is not 
constrained. 

In their phenomenological study, when comparing resummed and fixed-order NLO predictions, 
the authors of Ref. [52] find significant resummation effects. The reason is twofold. First, the cuts 
that are considered in Ref. ^2] are more restrictive, and thus resummation effects are made 
more relevant. Second, the comparison is done one order lower than ours (i.e. the NLL+NLO 
resummed prediction is compared to the fixed order NLO result) where we also find more significant 
distortions of the relevant kinematical distributions (see Figs. I MTOl) . 



4 Summary and outlook 

We have presented a calculation of the NNLL+NNLO cross section for Higgs boson production 
at the LHC, in the decay modes H — )■ 77, H — >■ WW — )■ Ivlv and H — )■ ZZ — )■ 4 leptons. The 
calculation takes into account some illustrative experimental cuts analogous to the ones designed 
to isolate the Higgs boson signal. 

Our calculation is implemented in the numerical program HRes [37j. The present version 
of the program includes the most relevant decay modes of the Higgs boson, namely, H — )■ 77, 
H — !■ WW — )■ Ivlv and H -^ ZZ — )■ 4 leptons. In the latter case it is possible to choose between 
H —7- ZZ —7- yU+/i^e^e~ and H — ?■ ZZ — )■ e^e~e^e~, which includes the appropriate interference 
contribution. The user can apply all the required cuts on the Higgs boson and its decay products 
and plot the corresponding distributions in the form of bin histograms. These features should 
make our program a useful tool for Higgs searches and studies at the Tevatron and the LHC. 

The calculations performed through HRes strictly implement the large Mt approximation. This 
is known to be a good approximation for the px spectrum of the Higgs boson, provided that px 
is not too large [pT^Mt) [55]. For very large transverse momenta the large- M^ approximation 
is bound to fail, since the QCD radiation accompanying the Higgs boson becomes sensitive to 
the heavy-quark loop. The inclusion of top and bottom mass effects up to 0{a^) in HRes is 
feasible and is left to future work. Another limitation of the calculation is that we completely 
neglect radiative corrections in the Higgs boson decay. The full QCD+EW corrections to the 
decay modes H — )■ WW{ZZ) — )■ 4 leptons are available ^6] and we plan to include these effects 
in a future version of the program. 
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